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AMI Galactic Plane Survey at 16 GHz: I - Observing, mapping and 
source extraction* 
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ABSTRACT 

The AMI Galactic Plane Survey (AMIGPS) is a large area survey of the outer Galactic plane to 
provide arcminute resolution images at milli-Jansky sensitivity in the centimetre-wave band. 
Here we present the first data release of the survey, consisting of 868 deg 2 of the Galactic 
plane, covering the area 76° g ( ^ 170° between latitudes of \b\ ^ 5°, at a central frequency 
of 15.75 GHz (1.9 cm). We describe in detail the drift scan observations which have been used 
to construct the maps, including the techniques used for observing, mapping and source ex- 
traction, and summarise the properties of the finalized datasets. These observations constitute 
the most sensitive Galactic plane survey of large extent at centimetre-wave frequencies greater 
than 1.4 GHz. 
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1 INTRODUCTION 

Large-area radio surveys contribute to our understanding of the 
Universe in numerous and diverse ways. Discoveries from these 
surveys have become key ingredients of moder n astrophysics : pul- 
sars, radio galaxies, quasars and more (see e.g. iLongaitll 19981) . For 
studies of our Galaxy, radio surveys are particularly beneficial as 
the longer wavelength radio emission does not suffer from the same 
extinction and opacity effects as optical and infra-red surveys and 
the dense regions of dust and gas which dominate the low-latitude 
Galactic plane become largely transparent, allowing us to study 
sources in these regions. However, the bulk of Galactic radio sur- 
veys are at frequencies at or below 1 .4 GHz and as such are nec- 
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essarily biased against objects whose spectra rise with frequency, 
such as dense star-forming regions. Two examples of the need for 
higher-frequency, centimetre-wave Galactic surveying are as fol- 
lows. 

The first is the hypercompact Hn (HCHu) region. Thought to 
indicate the earliest visible stage of massive star formation, these 
objects are two orders of magnitude more dense than the better 
known ultracompact (UCHu) region and have steeply rising spec- 
tra. HCHu regions were discovered serendipitously in observa- 
tions of UCHu, having been missed previously in their entirety by 
Galactic plane surveys concentrated at v < 5 GHz. The turnover 
frequency between the optically thick and thin regimes for ther- 
mal bremsstrahlung is a lin ear function of emission measure (e.g. 
iMezger & Henderson! 1 19671) causing such low frequency surveys 
(e.g. v < 5 GHz) to preferentially select against dense plasmas 
(« c < 10 11 m" 3 ). Such plasmas are not limited to HCHu regions 
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but also include a variety of other Galactic objects such a s massive 
stella r winds, ionised jets from young st ellar objects (e.g. lAngladal 
Il995h and young planetary nebulae (e.g. lBains et al ]|2009l) . 

The second is the anomalous microwave emission (AME), 
now being identified in an increasing number of Galactic ob- 
jects, that was missed in low fre quency Galactic su rveys. First 
identified by CMB experiments dLeitch et all 1 19971) as a large 
scale foreground contaminant this form of emission has since 
been demonstrated to exist i n more compact objects such as dark 
(e.g. ICasassus et al ] 120061 : lAMI Consortium: Scaife et al ] 120091: 



Scaife et al. 
Tibbs et all 



20101) and molecular clouds (Watson et al. I 120051 : 



20111) . Although multiple mechanisms have been pro- 



posed to explain A ME, dipole emission from rapidly rotating very 
small dust grains jDraine & Lazariar] Il998al . iDraine & Lazarianl 
Il998bl) is generally considered to be most likely. Such spinning 
dust emission has a sharply peaked SED with a maximum in the 
frequency range 10-50 GHz depending on grain size distributions. 

A current lack of surveys in this frequency range means 
that our knowledge of the overall properties of objects which ex- 
hibit emission from spinning dust, objects which are character- 
ized by dense plasmas, and indeed the global distribution of rising- 
spectrum emission in the Galaxy, is extremely poor. Those surveys 
which are available, s uch as the 9C Ryle Telescope s urvey (15 GHz; 
Waldram et al]|2003h . the GPA survey ( 14.35 GHz;lLangston et al.l 
200d) and the AT20G survey (20 GHz; iMiirphv et alj|20ich have 
provided us with tantalising insights into the high frequency Galac- 
tic plane, but there is a continuing need for higher sensitivity, reso- 
lution and sky area coverage at these frequencies. 

The interferometric Arcminute Microkelvin Imager (AMI) 
Galactic Plane Survey (AMIGPS) provides the most sensitive 
centimetre-wave Galactic plane survey of large extent at v > 
1.4 GHz. AMIGPS is a drift-scan survey of the northern Galac- 
tic plane at « 16 GHz, covering (in the first data release) the re- 
gion 76° < t < 170° and \b\ % 5°. The AMI Small Array (SA) 
has been used for the survey since its relatively large field of view 
(« 400arcmin 2 ) makes covering large areas feasible, and its short 
baselines mean that extended objects, very common in the Galaxy, 
are at least partially observable. The resolution of the survey is 
ss 3 arcmin and the noise level is a 3 mjy beam 1 away from bright 
sources. 

This paper focuses on the techniques employed for observing 
(Section[3]l, mapping (Section[4]l and source extraction (Section[5j 
in the AMIGPS. The positional and flux density calibration accu- 
racy of the survey are also tested in Section[6] and in Section[7]the 
maps and catalogue are described. In a following paper, hereafter 
Paper II, the first results from the survey, including the follow-up 
of rising-spectrum objects in order to detect U/HCH11 regions, will 
be presented. 



will be used as calibration accuracy checks in Section|6]and will be 
described fully in Paper II. 

The SA is an interferometer array comprising ten 3.7-metre- 
diameter, equatorially-mounted dishes, with a range of baselines 
of 5 - 20 m. It operates at frequencies over 13.9 - 18.2 GHz with 
the passband divided into six channels of 0.72-GHz bandwidth. 
It has a primary beam at the central frequency of 15.75 GHz of 
as 20 arcmin full width at half-maximum (FWHM) and a typical 
synthesised beam FWHM of a 3 arcmin (this varies depending on 
the precise «v-coverage of any observation). The telescope mea- 
sures a single, linear polarisation (Stokes 1 + Q) and has a flux- 
density sensitivity of ss 30 mjy s~' /2 . It is sensitive to angular scales 
of up to a 10 arcmin (depending on the Mv-coverage). 



3 OBSERVATIONS 

The AMIGPS is observed in drift-scan mode, in which the SA is 
pointed at a fixed azimuth and elevation and observes a narrow 
strip as the sky drifts past. In practice, to enable re-observation of 
strips as necessary, the telescope is driven very slowly to maintain 
a constant J2000 declination. In order to perform phase calibration, 
bright nearby point sources sel ected from the Very Long Baseline 
Array Calibrator Survey (VCS. lBeaslev et al Eool) were observed 
for 400 seconds at 30-minute intervals during each scan. Strips are 
observed at a separation of 12 arcmin in 6, corresponding to the 
35% point of the power primary beam, i.e. at distance x from the 
centre where exp(-x 2 /(2cr 2 )) = 0.35, assuming the beam is Gaus- 
sian with width tr. This produces a very even noise level across the 
combined map, with a variation of a 3% between the centre of a 
declination strip and the point halfway between declination strips. 
The noise level in the survey is « 3mJybeairr' away from bright 
sources and is as low as a 1 mjy beam 1 at some points. 

This data release consists of observations above 6 = 40° and 
between b « ±5°; a later data release will extend the coverage to 
S > 20°, corresponding to 53° g I g 76° and 170° g t g 193° . 
The coverages of some other, currently available Galactic plane 
surveys along with their resolutions and noise levels are shown 
in Table [T] and some of these are illustrated in comparison to the 
(full) AMIGPS in Fig. □ The AMIGPS is the first survey at this 
frequency to achieve similar coverage area, resolution and noise 
level to lower freque ncy surveys such as the Canadian Galactic 
Plane Survey (CGPS; iTavlor et alj|2003l) ; earlier surveys have ei- 
ther been wide and shallow with lower resolution (e.g. the GPA), 
or narrower, with comparable resolution but still more shallo w than 
the AMIGPS (e.g. Nobeyama at 10GHz; lHaiida et al.lll987l) . 

The observations for the first data release were performed be- 
tween 22 Jun 2010 and 4 Nov 2011. Approximately two thirds of 
the strips were observed multiple times in order to improve the 
noise level, resulting in a total observing time of « 1200 hours. 



2 THE ARCMINUTE MICROKELVIN IMAGER SMALL 
ARRAY 

The AMI dAMI Consortium: Zwart et al]|2008h was designed as a 
Sunyaev-Zel'dovich (SZ) effect instrument, and consists of two 
separate instruments: the Small Array (SA), optimised for observ- 
ing extended SZ decrements on arcminute scales, and the Large 
Array (LA) with higher resolution (30arcsec), used for charac- 
terising and subtracting point source foregrounds from SA data. 
The AMIGPS was carried out solely with the SA, although some 
follow-up observations were also made using the LA; these latter 



4 DATA REDUCTION AND MAPPING 

Data reduction was performed using the local software tool re- 
duce, which flags interference, shadowing and hardware errors, 
applies phase and amplitude calibrations, and Fourier transforms 
the lag correlator data to synthesize the frequency channels, be- 
fore outputting to disc in mv-fits format. Flux calibration was per- 
formed using short observations of 3C48, 3C286 or 3C147 near 
the beginning and end of each run. The assumed flux densities for 
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Figure 1. Coverage of the full AMIGPS compared to other northern Galactic plane surveys of similar area and/or resolution and noise level. The AMIGPS 
boundaries are shown as a solid black line, CGPS (408, 1420 MHz) as a dashed black line, Effelsberg (1.4, 2.7 GHz) as a dotted black line, 7C(G) (151 MHz) 
as a solid cyan line, Nobeyama (10 GHz) as a dashed cyan line, and CORNISH (5 GHz) as a dotted cyan line. 



3C286 were calculated from Very Large Array total-intensity mea- 
surements provided by R. Perley (priv ate communication), and are 
consistent with the iRudv et al.l d 19871) model of Mars transferred 
onto absolute scale, using results from the Wilkinson Microwave 
Anisotropy Probe satellite. The assumed flux densities for 3C48 
and 3C147 are based on long-term monitoring with the AMI SA 
using 3C286 for flux calibration (see Table [2j. A correction for 
each antenna for changing weather and radiometer performance is 
also applied using a noise-injection system, the 'rain gauge' (see 
AMI Consortium: Zwart et al ]|2008l for more details). 

Many of the automatic flagging routines used to excise inter- 
ference from AMI data rely on the amplitude of the astronomical 
source being observed remaining constant throughout the obser- 
vation. This is not the case for drift-scan data, where sources drift 
through the primary beam. It was found that the automated routines 
(designed in particular to remove interference spikes) were partially 
flagging out bright sources, so that their flux densities measured 
from the final maps were systematically lower than expected com- 
pared to flux densities from tracked observations. To overcome this, 
an iterative scheme was introduced in which the data are reduced as 
usual, then the visibility data are averaged over channels and base- 
lines and searched for remaining peaks, indicating the presence of 
a bright source. The data are then re-reduced, with the time ranges 
within which bright sources are found being excluded from the in- 
terference flagging routines. 



In order to process the continuous drift-scans, pointing cen- 
tres were set up spaced lOarcmin apart in RA and each sample 
was phase-rotated to the closest pointing centre. The SA integration 
time is one second, so this results in 600/(15 cos(<5)) samples per 
pointing centre, which yields approximately 70 samples at 6 = 55°. 
The data were then exported in multi-source mv-fits format. The 
in-house software package Fuse was used to concatenate visibil- 
ity data from different observations of the same declination strip, 
weighting each observation according to its noise level; variability 
of sources has not been considered (see Section. [6~2t . 

4.1 Mapping 

The pointings were imaged separately in Aipflwith 128x128 pixel 
maps, where the pixels are 20x20 arcsec 2 in size. Natural weighting 
was used to maximise signal-to-noise, and all six frequency chan- 
nels were imaged using a multi-frequency synthesis; as a conse- 
quence of different flagging of the channels, the effective frequency 
will vary slightly between pointings. Individual channel maps were 
not produced. An automated CLEANing process was used to place 
a 6 x 6 pixel Clean box around the brightest pixel if it had a flux 
density greater than 200 mjy, and the images were CLEANed to three 

1 astronomical image processing system - www.aips.nrao.edu/ 
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Table 1. Coverage, resolution and noise levels of selected Galactic plane 
surveys. The noise level marked with (*) is actually a detection limit. 



Telescope/ 


Frequency 


Coverage 


Resolution 


Noise level 


Survey name 


(GHz) 


(deg 2 ) 


(arcmin) 


(mjy beam - ') 


7C(G)" 


0.151 


1700 


1.17 cosec((5) 


40 


Al 20G" 


20 


20086 


1 .7 


10 


CGPS C 


1.42 
0.408 


1500 


1 cosec(<5) 
3.4cosec(<S) 


0.23 
3 


CORNISH f/ 


5 


1 10 


0.017 


0.4 




4.875 1, 


125 


2.6 


120(*) 


cttelsberg 


1 Af'S 

n nhj 
Z. / 


2400 


Q ^ 
J. j 


JU 


GPA J 


8.35 
14.35 


2700 


10 

7 


230 
800 


MAGPIS* 


1.42 


43.2 


0.083 


0.2 


Nobeyama' 


J u 


1 OJ 




jj 




0.96 




4 x 75 


60 


RATAN"' 


3.9 


400 


1 x39 


10 




11.2 




0.35 x 14 


100 


Stockert" 


2.72 


10200 


18 


140 


VGPS" 


1.42 


< 200 


1 


2 


VLA 


5'' 


40 


0.07 


2.5-10 


1.42« 


224 


0.07 


10 


VSA r 


30 


152 


13 


90 



" Vessev& Green (1998); 6 Murph y et al. (2010) ; c Taylor et al. (2003) : 
" Purcell. Hoare. & Dia mond (2008); e Alten hoff et al. (1979) ; 
f Reich. Reich. & Fiirst ( 1990) ; s Reich. Reich. & Fiirst ( 1997) ; 
Reich etal. 0984) : ' Reich et al. 0990) : 1 Langston et al. (2000) : 
* Helfandet al. (2006): ' Handa et al. (1987): "' Trushkin (1998): 



" Rett et al.TI987 ): " Shi et al (20U6): Becker et al. (1994); 
^ zoon ematkerman i et at. ( 1^90) : ' Todorovic et al. 



Table 2. Assumed I + Q flux densities of 3C286, 3C48 and 3C147. 



Channel 


v/GHz 


s 3C286 /Jy 


5 3C48 /Jy 


S 3C147 /Jy 


3 


13.88 


3.74 


1.89 


2.72 


4 


14.63 


3.60 


1.78 


2.58 


5 


15.38 


3.47 


1.68 


2.45 


6 


16.13 


3.35 


1.60 


2.34 


7 


16.88 


3.24 


1.52 


2.23 


8 


17.63 


3.14 


1.45 


2.13 



times the r.m.s. noise on the dirty image; the box was removed, and 
the CLEANing was continued to the same flux density level. Each 
component map was CLEANed using an elliptical Gaussian fitted to 
the central region of the dirty beam as the restoring beam. As a re- 
sult, the restoring beam for each component map is slightly differ- 
ent. The distribution of synthesised beam sizes is shown in Fig. [2] 
The noise on each component map was estimated using the 
Imean task over the whole map, which fits a Gaussian centred on 
zero to the distribution of pixel values, discarding outliers. Fig. [3] 
shows a typical beam-corrected pointing map at S ~ 55° with its 
restoring beam, along with the nv-coverage for the pointing centre. 



4.2 Primary-beam correction 

The S A primary beam is usually approximated by a Gaussian fitted 
to the central lobe of the beam, with a FWHM of 19.6 arcmin at the 



1 



| | Major axis 
Minor axis 



150 200 250 300 
Synthesised beam FWHM (arcscc) 



Figure 2. The distribution of synthesised beam major and axis FWHMs for 
the pointing centres which make up the drift-scan maps. 



central frequency. However, primary beam correction of drift-scan 
maps is complicated due to the continuous nature of the scan: each 
pointing does not consist solely of data taken towards its centre, 
but rather of data taken towards a series of points along its RA axis, 
corresponding to data taken at different times. The primary beam 
was therefore calculated as a weighted average of beams centred at 
each of the constituent points, i.e. for any pixel in the map 



primary beam 



Zf =1 w,exp(-i^) 



(1) 



where N is the number of samples constituting the pointing, w, = 
l/cfmsi is the weight of the f'th sample (i.e. the sum of weights 
for all baselines and all channels contributing to a one-second sam- 
ple) where cr rms- i is the r.m.s. noise on the sample, 2cr V(21n(2)) 
is the FWHM of the SA primary beam (19.6 arcmin), and A; = 
y/(x - x, ) 2 + y 2 is the separation of the pixel from the pointing cen- 
tre of the sample, where (x, y) is the pixel location and (x,, 0) is the 
pointing centre of the sample along the RA axis. The pixel value 
is then divided by the weighted-average beam for that pixel; pixels 
with a weighted- average beam of < 0.1 are blanked. This has the 
effect of elongating the beam along the RA axis compared to the 
normal SA primary beam correction; this is illustrated in Fig. [4] 

This beam correction was applied to each of the individual 
pointing maps outputted from Aips and a noise map for each point- 
ing was also produced, which is the inverse of the primary beam 
scaled by the r.m.s. noise value of the map outputted from Imean. 



4.3 Map combination 

Finally, the individual beam-corrected pointing maps were added 
together, weighting each pixel by the inverse of its variance calcu- 
lated from the noise ma p, into larger continuo us maps using the in- 
house software Profile jGrainge et al. 2002). Corresponding con- 
tinuous noise maps for use in source-finding were also produced 
in the same way from the noise maps for the individual pointing 
centres; these are found to provide an accurate representation of 
the noise, except around bright sources as discussed in Section PJSl 
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Figure 3. The beam-corrected map (top) and wv-coverage (bottom) for a 
typical drift-scan pointing centre at 8 ss 55° . The la" map noise at the centre 
of the map is 2.72mjybeam~', as estimated by Imean; the noise increases 
away from the centre (up to a factor of ten). The map has contours at ±2— 
10o"; solid contours are positive and dashed contours are negative. The syn- 
thesised beam is shown in the bottom left-hand corner. The grey-scale is in 
units of mJy beam -1 . 
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Figure 4. An example of the drift-scan primary beam. This is the weighted 
average of the SA primary beam centred at each of the samples taken along 
the RA axis. The figure on the top shows the weighted-average primary 
beam for the same pointing shown in Fig. [5] the bottom figure shows the 
weighted-average primary beam divided by the usual SA primary beam for 
a typical, non-drift-scan pointing. The coordinates of the samples constitut- 
ing the pointing are overplotted in black. The diameter in RA of the pointing 
is « 37 arcmin, whereas the diameter of a typical, non-drift-scan pointing 
cut off at the 10% point is x 35 arcmin. 



Fig.[5]shows an example noise map section illustrating the variation 
in noise level across a typical map. All maps were also regridded 
into Galactic co-ordinates using the Aips task regrd. 



5 SOURCE EXTRACTION 

AMI Consortium: Franzen et al. | (Em ]]) describes the source ex- 
traction methodology used for the 10C survey; a similar process 
was used for the AMIGPS. Source finding was carried out over the 
combined maps to search for peaks greater than 5cr, where cr is the 
r.m.s. noise value read from the corresponding pixel on the noise 
map. A peak position and flux density value is measured by inter- 
polating between the grid points. An initial estimate of the inte- 
grated flux density and source size is also calculated by integrating 
contiguous pixels down to 2.5x the local thermal noise level, and 
sources are identified as overlapping if the integration area contains 



more than one peak > 5cr. This information is used to fit an ellip- 
tical Gaussian to each source in an automated fashion, using the 
Aips task Jmfit. Overlapping sources are fitted simultaneously. As 
Jmfit can only fit up to four sources simultaneously, in a limited 
number of cases where more than four sources were identified as 
overlapping by the source-finding software sources were regrouped 
manually into groups of four or less. 



5.1 Source size and classification 

The deconvolved source size is calculated by Jmfit using the syn- 
thesised beam size at the pointing with the highest weight at the 
position of the source. This size is used to classify the source as 
point-like or extended to the SA beam, followin g the method de- 
scribed in lAMI Consortium: Franzen et al.lf201ll) . scaled to the SA 
beam size. A source is classified as extended if the fitted major axis 
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Figure 5. A typical noise map illustrating the variation in noise level across 
the map. The grey-scale is in mJy beam - 1 and is truncated at 6 mJy beanT 1 
to show the low-level variation; the highest noise level in the area shown is 
« 10mjy beairT 1 in the north-eastern corner. Crosses mark the positions of 
sources with peak flux densities > 50 mJy beam - 1 , around which it can be 
seen that the noise level increases. Away from the bright sources, the noise 
level is g 3 mJy beam -1 . 
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Figure 6. Pixel values (solid black lines) interpolated through map points 
in lines intersecting the bright, central source in Fig. [JJ and the spuri- 
ous sources around it; the mean noise- and 5cr-detection-levels (solid and 
dashed cyan lines) from the noise map; the fall-off law and exclusion 
zone radius for this source (cyan curved and vertical dotted lines); and the 
S peak.bright/10 cutoff line (dashed black line). It can be seen that the noise 
outside the exclusion zones is well represented by the map noise, however 
closer to the central source the noise is elevated and the 5<x cutoff is not 
high enough. The conservative Sp e ak,bright/10 cutoff excludes the spurious 
detections. 



size e ma j > e crit , where 

= f 3.02> maj p- 11 - if 3.0£ maj p- 1 ' 2 > 100.0 arcsec, 
ecnt \ 100.0 arcsec otherwise, K) 

where p is the signal-to-noise ratio and £> raa j is the synthesised-beam 
major-axis size. 

If a source is classified as extended, its integrated flux density 
outputted from Jmfit is considered to best represent its total flux 
density; otherwise the peak flux density is considered to provide a 
more accurate measurement. 

The error inherent in using the beam from the pointing with 
the highest weight at the position of the source for source extraction 
has been investigated by remapping a section of the survey with 
identical restoring beams for all pointing centres. The flux densities 
derived from this map were compared with the catalogue values 
for sources which lie between pointing centres with different beam 
shapes and sizes. For point-like sources, the difference in the flux 
density is g 1 % and is considered to be negligible. For extended 
sources it is % 5%, so a conservative extra 5% error on the flux 
density is added in quadrature (see Section[6j2}. 

Although a Gaussian is a reasonable approximation to the 
shape of many sources, clearly in the Galactic plane there are many 
sources which are not Gaussian in shape. Integrated flux densities 
should therefore be used with caution. The '% 2 ' statistic (actually 
the sum of squared residuals) outputted from Jmfit is included in 
the catalogue as an indication of the goodness of fit. It should be 
noted that integrated flux densities will also be affected by flux loss 
when the source is extended on scales that are large compared to 
the SA synthesised beam. 

5.2 Spurious source exclusion 

An implausibly large number of sources are frequently detected 
in the vicinity of bright sources - these are likely spurious and 



are caused by residual amplitude and phase errors in the data and 
unCLEANed sidelobes. In order to prevent these from contaminating 
the catalogue, 'exclusion zones' were applied to sources with peak 
flux density > 50mJy beam -1 . The radii r E of the exclusion zones 
are determined by the peak flux density S P eak,bright of the source as 

/ \ 1/3 

r E = 18 (SpeaiLbrjeht/SOOmJy) arcmin. This was chosen empiri- 
cally to describe the fall-off in the elevated, non-Gaussian noise 
around bright sources, illustrated in Fig. [6] Within the exclusion 
zones, only 'sources' with peak flux density 5 pca j; > Speak.bright/10 
are retained. The factor of ten was conservatively chosen by eye 
to retain most of the sources which appear to be real, while ex- 
cluding as many spurious sources as possible. There may be some 
real sources which are excluded by this procedure; the implications 
of this for the completeness of the catalogue will be discussed in 
Paper II. Fig. [7] illustrates the exclusion zones around two bright 
sources. 



6 CALIBRATION ACCURACY CHECKS 

As many bright sources found in the Galactic plane are well known 
and frequently used as phase-calibrator sources by AMI, it is pos- 
sible to use them to check both the positional and flux calibration 
accuracy of the drift-scan survey. A follow-up campaign with the 
AMI LA has also been conducted on so urces identified as h aving 
rising spectra with respect to the NVSS ^Condon et alJI 199*3) (this 
will be described in detail in Paper II). Since the positional accu- 
racy of the LA is much greater than that of the SA (the synthesised 
beam is «30 arcsec), these pointed follow-up observations provide 
an additional check on the calibration accuracy. 
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22 54 52 50 48 

RIGHT ASCENSION (J2000) 

Figure 7. A section of the map illustrating the spurious source exclu- 
sion method. The grey-scale is in mJybeairr 1 and is truncated to show 
the fainter sources; the flux densities of the brightest and second brightest 
sources are «1700mJy and 50mJy respectively. The contour levels are be- 
tween ± 1 00 mJy beairT 1 in steps of 1 mJy beam" 1 (it is not possible to use 
<x contours since the noise level varies across the map); solid contours are 
positive and dashed contours are negative. Exclusion zones are shown as 
circles around the bright sources. Source detections are marked by X, and 
'sources' detected but excluded by +. The synthesised beam at the position 
of the brightest source is shown in the bottom left hand comer. 



Table 3. Mean RA and S position offsets for high SNR sources in the drift- 
scan catalogue (all in arcsec). Consistency is checked by using the offsets 
from VCS catalogue and LA positions separately, and combined, and by 
changing the minimum SNR. 



SNR 


Offset 


Number of 


Mean RA 


Mean S 


limit 


from 


sources 


offset 


offset 




VCS 


56 


0.9 ± 0.5 


-0.4 ± 0.3 


50 


LA 


18 


-1.1 ±0.7 


-0.9 ± 0.7 




Combined 


74 


0.5 ± 0.4 


-0.5 ±0.3 




VCS 


30 


0.6 ±0.7 


-0.6 ± 0.4 


100 


LA 


5 


-0.5 ± 1.5 


-1.0 ± 1.0 




Combined 


35 


0.5 ± 0.4 


-0.65 ± 0.3 




2 
RA offset / a RA 



6.1 Positional accuracy 

6.1.1 Point-like sources 

The catalogue of source positions (for point-like sources only) de- 
rived from the survey maps was matched t o the milliarcsecond - 
accurate positions from the VCS catalogue dBeaslev et alj|2002r) . 
resulting in 125 matches with signal-to-noise ratio (SNR) in the 
drift-scan maps ranging from * 8 to 600. In addition, the posi- 
tions derived from follow-up observations of objects that were also 
point-like to the LA were compared to the drift-scan catalogue po- 
sitions, resulting in 270 additional matches (not matched to a VCS 
source) with SNR in the drift-scan maps ranging from a 5 to 400. 

The errors cr^ and cr B in RA and 8 for a point source are 
assumed to be given by 



' RA 

2 



4a + °m sin 2 (0) + of n cos 2 (0) 

2 



(3a) 

' e ~ 4 + °m cos-(0) + cr^ sin z (0), (3b) 

where egAort are the r.m.s. calibration errors in RA and 6, o"Morm 
are the noiselike uncertainties parallel to the synthesised beam ma- 
jor (M) and minor (m) axes, and <p is the position angle of the beam. 
We assume the noiselike uncertainties are given by 

c ™°' m = vwsnr' (4) 

where 0M Orra are the major and minor FWHM of the synthesised 
beam. 

To test for systematic RA and 5 offsets, the mean offsets be- 
tween both the AMIGPS and VCS catalogue and AMIGPS and LA 



Figure 8. RA and S offsets normalised by their calculated errors for all 
sources which are matched to a VCS source with well-known position 
(black crosses) or have been followed up with the LA (cyan dots). The esti- 
mated 1 and 3<r error circles are also shown. 



positions were calculated separately and as a single group, and by 
selecting sources with SNR >50 and SNR >100 in the drift-scan 
maps. These are listed in Table [3] and are all consistent with zero 
within < 2.5cr, so we assume no systematic offset in RA or 5. 

To determine the r.m.s. calibration errors, e RA and e 5 were 
varied until 99.7% of the sources with VCS positions had offsets 
within 3cr calculated from Equation [3] This gave era = 2.6 arcsec 
and e 6 = 1.7 arcsec. Fig. [8] shows the positional offsets for all 
sources in both datasets, normalised by the calculated error. They 
agree well, with 99% of all offsets lying within the 3cr circle (the 
statistics are only approximately Gaussian). 



6.1.2 Extended sources 

For sources that are extended given the SA beam, the positional 
uncertainty is calculated slightly differently. The errors in RA and 
6 are given by 

°~RA ~ e RA + °7,RA P 3 ) 

o-] = 4 + o-], 5 , (5b) 
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where the o- JRAoT]tS terms are the errors outputted by the Aips fit- 
ting task Jmfit, which folds in an estimate of the noiselike error as 
well as the error associated with the fit. 

6.2 Flux calibration accuracy 

We assume flux calibration errors are given by 

Og ^ = (0.05 Speak) 2 + o" 2 for a point-like source (6a) 

crj; = 2 (0.05 5 j nt ) 2 + a 2 for an extended source, (6b) 

where S pea ]t is peak flux density and 5 j nt is integrated flux density. 
This error estimation comprises a 5% calibration uncertainty (in- 
cluding rain-gauge correction) and a noise-like error cr which for a 
point-like source is the r.m.s map noise measured from the CLEANed 
map, and for an extended source is the error outputted from Jmfit 
which also folds in an estimation of the fitting error. The error for 
an extended source also contains an extra 5% error due to the un- 
certainty in the beam shape, as discussed in Section|5] 

At 16 GHz, intrinsic source variability is important. 
lAMI Consortium: Franzen et al.l d2009h find that of 93 extra- 
galactic sources monitored with the AMI SA for periods between 
one and 18 months, «50% are variable above the flux density cali- 
bration uncertainties and 15% are variable at a level of more than 
20%. Variability must therefore be considered when attempting to 
test the flux calibration accuracy. 

6.2.1 NGC7027 

NGC 7027 is a planetary nebula lying within the drift-scan sur- 
vey area_and_Jo^_rjres^nt_rjur^s£sessentially non- variable (see, 
e.g. IZiilstra. van Hoof. & Pe rlev 2008). It is also frequently mon- 
itored by AMI with tracked observations so an accurate flux den- 
sity at 16 GHz can be calculated for comparison. Using data taken 
between 2007 and July 2012 with the SA, the 16 GHz flux of 
NGC 7027 is 5 .4 Jy. The drift-scan flux for NGC 7027 is 5. 1 ±0.3 Jy, 
agreeing with the tracked value to within 1.2cr. 

6.2.2 3C48 

3C48 is one of the primary calibration sources used by AMI and 
is also known to be non-variable. Separate drift-scan observations 
were made of an area around it between Mar and Dec 2010 as an 
initial test of the drift-scan pipeline. These observations were re- 
duced both in the standard pipeline, which uses the closest primary 
calibrator observations in time including 3C48, as well as using 
only 3C286 as a primary calibrator. 3C48 has a flux density of 
1.64Jy at 16 GHz; the drift-scan flux density is 1.60±0.08 Jy, us- 
ing 3C48 and 3C286 as primary calibrators, and 1.63±0.08 Jy us- 
ing only 3C286. Both values are consistent with each other and are 
within 0.5cr of the nominal value. 

6.2.3 Concurrent observations 

Since AMI is continually observing phase calibrators for many of 
its observations, there is a high probability of quasi-simultaneous 
tracked observations existing of bright compact sources, mostly 
extra-gal actic, seen in the drift-scan survey. E xtrapolating from 
Fig. 3 of lAMI Consortium: Franzen et alj d2009t) which shows the 
variability index for extra-galactic sources at 15 GHz as a function 
of time, 10 days was chosen as an interval within which source 




100 1000 10000 



Mean quasi simultaneous tracked flux density (mJy) 

Figure 9. Drift-scan flux densities compared to the mean flux from tracked 
SA archival observations within 10 days. The black solid and dotted lines 
show a one-to-one correspondence and ±5% flux-calibration uncertainty. 



variability should be small. Since the drift-scan survey also con- 
sists of multiple observations at different dates, each observation 
which contained a potential match within +10 days was reimaged 
separately and source-finding was done on the individual declina- 
tion strips. Any archival SA tracked observations within ±10 days 
of drift observations of matching sources were averaged and com- 
pared with the individual drift-scan values. Fig. [9] illustrates the 
comparison between the peak flux densities of these sources; 93% 
of the drift-scan flux densities are within 3tr of the mean archival 
flux. 

The three outliers had lower drift-scan flux densities than the 
mean archival flux density and were found to lie near the edge of 
the declination strip, where phase errors are expected to have the 
greatest effect. In each case, the source appears near the centre of 
the adjacent strip, which was observed a day later. When creating 
the final combined map, the pixels nearer the centre of individual 
pointings are given greater weight, so the discrepant flux densities 
will be down-weighted. The flux densities for these sources derived 
from raster maps produced from observations close in time agree 
with the mean archival flux to within l<x. 

It is common for survey flux densities to be slightly suppressed 
due to small phase errors shifting the positions of sources which lie 
away from the pointing centres in th e constituent maps (see, e.g. 
I AMI Consortium: Davies et alj201 lh . The concurrent observations 
were tested for this effect, but the median percentage difference 

((Smean.tracked ~ S drift) /Sm«n,tracked) was found t0 be onl y ~ 2% > so 

the AMIGPS flux densities were not adjusted for this effect. 
6.2.4 Non-concurrent observations 

A final check of the flux calibration accuracy can be made by 
comparing the LA follow-up flux densities to the drift-scan flux 
densities for sources that are found to be point-like to the LA, 
although these observations are widely spaced in time (by up to 
as 1.5 years). Very little is known about variability statistics in the 
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Figure 10. Drift-scan flux densities compared to the LA follow-up flux. The 
black solid line shows a one-to-one correspondence. 



Galactic plane at cm-wavelength. However, some idea of the ex- 
pected number of variable sources can be obtained using results 
from the 5-GHz Galactic plane variability study by l lBecker et al.l 
l201Ch . where a 8% of sources detected in the flux density range 
from 1-100 mJy between b ss ±1.0° were found to be variable on 
a time scale of years or shorter when no correction for the inclusion 
of the extra-galactic source population was applied. 

Fig. [10] shows the comparison between the LA and drift-scan 
survey peak flux densities. 87% are within 3cr, taking into account 
the LA errors which are generally smaller than the drift-scan errors 
and are not plotted for clarity. The remaining 13% seems consis- 
tent with the 8% of sources predicted to be variable, given that no 
correction for differences in frequency, flux density range, Galac- 
tic latitude distribution or bias due to selecting for rising spectrum 
sources has been attempted. The apparent bias towards higher drift- 
scan flux densities at the lower end of the flux density scale is likely 
an Eddington bias caused by low-SNR sources selected from the 
AMIGPS map being more likely to occur on positive noise bumps. 



7 DATA PRODUCTS 
7.1 Raster maps 

The field is divided into 38 square maps of side 6°, and are given 
names constructed from the Galactic coordinates of their centres, 
eg G78.0-2.2. These are shown in Fig.QT] The centres are spaced 
by 5° in longitude, and 4.4° in latitude, and start at I = 78.0°, b = 
-2.2°. 

These raster maps, along with the corresponding noise maps, 
are available from http://www.mrao.cam.ac.uk/surveys/AMIGPS/ 
Also included with the maps as an extra fits data-cube are tables 
giving the synthesised beam sizes appropriate to each pixel (i.e. the 
synthesised beam belonging to the pointing with the highest weight 
at that pixel). 

Fig.[J_2]shows an example 6 deg 2 map, with annotations mark- 
ing the sources detected within it. Also shown for comparison is 
a CGPS total intensity 1.4 GHz map showing the same region. It 
can be seen that many sources detected by CGPS are also detected 
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Tycho 



I L_ I I I 

115 116 117 118 119 120 121 

GALACTIC LONG. 

Figure 12. An example AMIGPS raster map centred at I = 1 18.0°, b = 2.2° 
(top) and the corresponding CGPS total intensity map at 1.4 GHz (bottom). 
The AMIGPS map has source detections marked with x; the grey-scale is 
in truly beam -1 and is truncated at -20 and +50mJy beam -1 to show the 
fainter features. The grey-scale of the CGPS map is in K and is truncated at 
15 K for the same reason. Some well-know n supernova remnants (SNR) a nd 
Hn regions visible in the map are labelled jGreenl2009USharplesslll959l) . It 
can be seen that the AMIGPS sees many features common to the CGPS, 
however the larger-scale features such as the SNR Gl 16.5+1 . 1 are resolved 
out. 



by the AMIGPS; however some larger-scale features such as the 
supernova remnant Gl 16.5+1.1 are resolved out. 



7.2 Source catalogue 

An abbreviated sample of the catalogue containing the first ten 
sources detected in Fig. [12] is shown in Table [4] The com- 
plete source list, which contains 3503 entries, is available from 
http://www.mrao.cam.ac.uk/surveys/AMIGPS/ For each source, 
the catalogue contains: 

• A source name, constructed from the 12000 RA and 6 coordi- 
nates of the source. 
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Galactic longitude (degrees) 

Figure 11. The positions of the raster maps in Galactic coordinates. The solid black line marks the extent of the data, the dotted and dashed lines show the 
boundaries of the raster maps and the crosses mark the centres of the maps. 



• The peak RA, 6, flux density and associated errors (these are 
the appropriate quantities to use for point-like sources). 

• The fitted centroid RA and 6, integrated flux density and as- 
sociated errors (these are the appropriate quantities to use for ex- 
tended sources). 

• The critical source size as defined in Eqn. [2] and the decon- 
volved source major and minor axis sizes and position angle. A 
deconvolved size of 0.0 indicates that the source was not found to 
be wider than the synthesised beam in the major or minor axis di- 
rection. 

• The sum of squared residuals from the fit. 

• The source classification (point-like or extended). 



8 CONCLUSIONS 

The Galactic plane between b ss ±5° has been surveyed us- 
ing the interferometric AMI SA at « 16 GHz, to a noise level of 
~ 3 mjy beam -1 at a 3 arcmin resolution. This is the most sensitive 
and highest resolution Galactic plane survey at cm-wave frequen- 
cies above 1.4 GHz. 

(i) 868 deg 2 of the Galactic plane have been surveyed and a cat- 
alogue of 3503 sources produced. This is the first data release of 
the AMIGPS. 

(ii) As part of creating the AMIGPS, we have developed an au- 
tomated pipeline to produce maps from data taken in drift-scan 
mode, accounting for the presence of bright sources. 

(iii) The source extraction techniques developed for the 10C sur- 
vey have been applied to maps at different resolution and regions 
of the sky with many extended sources present. 

(iv) In testing the flux calibration of the survey by comparing 
source flux densities derived from the AMIGPS to tracked obser- 
vations of both extra-galactic and Galactic sources taken with the 
AMI SA and AMI LA, we find that the AMIGPS flux calibration is 
accurate to within 5%. 

(v) The r.m.s. positional accuracy of the survey, assessed by 
comparing positions derived from the AMIGPS with well-known 



source positions from the VLBA calibrator survey and with AMI 
LA follow-up positions, is 2.6arcsec in RA and 1.7 arcsec in S. 

In a following paper the first results from the survey will be 
presented, and in a future data release the survey will be extended 
to6> 20°. 
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